A constrained polynomial regression procedure for estimating the local False Discovery Rate

<p>Abstract</p> <p>Background</p> <p>In the context of genomic association studies, for which a large number of statistical tests are performed simultaneously, the local False Discovery Rate (<it>lFDR</it>), which quantifies the evidence of a specific gene association with a clinical or biological variable of interest, is a relevant criterion for taking into account the multiple testing problem. The <it>lFDR </it>not only allows an inference to be made for each gene through its specific value, but also an estimate of Benjamini-Hochberg's False Discovery Rate (<it>FDR</it>) for subsets of genes.</p> <p>Results</p> <p>In the framework of estimating procedures without any distributional assumption under the alternative hypothesis, a new and efficient procedure for estimating the <it>lFDR </it>is described. The results of a simulation study indicated good performances for the proposed estimator in comparison to four published ones. The five different procedures were applied to real datasets.</p> <p>Conclusion</p> <p>A novel and efficient procedure for estimating <it>lFDR </it>was developed and evaluated.</p

Long-term in vivo imaging and measurement of dendritic shrinkage of retinal ganglion cells

PURPOSE. To monitor and measure dendritic shrinkage of retinal ganglion cells (RGCs) in a strain of transgenic mice (Thy-1 YFP) that expresses yellow fluorescent proteins in neurons under the control of a Thy-1 promoter. METHODS. A total of 125 RGCs from 16 eyes of Thy-1 YFP transgenic mice were serially imaged with a confocal scanning laser ophthalmoscope for 6 months after optic nerve crush. Quantitative analysis of cell body area, axon diameter, dendritic field, number of terminal branches, total dendritic branch length, branching complexity, symmetry, and distance from the optic disc was used to characterize the morphology of RGCs, describe the patterns of axonal and dendritic degeneration, identify the morphologic predictors for cell survival, and estimate the rate of dendritic shrinkage. RESULTS. RGC damage was observed prospectively to begin with progressive dendritic shrinkage, followed by loss of the axon and the cell body. In a small proportion of RGCs, progressive axonal changes including fragmentation, beading, retraction, and bulb formation were also observed. RGCs with a larger dendritic field and a longer total dendritic branch length in general have a better survival probability. The rate of dendritic shrinkage was variable with a slower rate observed in cells having a larger dendritic field, a longer total dendritic branch length, and a greater distance from the optic disc. CONCLUSIONS. Estimating the probability of RGC survival and measuring the rate of dendritic shrinkage could become a new paradigm for investigating neuronal degeneration and evaluating the response of neuroprotective treatment. © 2011 The Association for Research in Vision and Ophthalmology, Inc.postprin

The role of diurnal cycle in subduction/obduction

Author Posting. © The Author(s), 2010. This is the author's version of the work. It is posted here by permission of Springer for personal use, not for redistribution. The definitive version was published in Journal of Oceanography 67 (2011): 273-279, doi:10.1007/s10872-011-0025-4.The annual subduction/obduction rate can be calculated in Lagrangian and Eulerian coordinates. In previous studies such calculations were primarily focused on the case with the seasonal cycle only. By extending these calculations to the case including the diurnal cycle of mixed layer depth, the annual subduction/obduction rate can be greatly increased.LLL and FW were supported by National Natural Science Foundation of China under Grant 40906007 and 40890150

Spatial constancy of attention across eye movements is mediated by the presence of visual objects

Recent studies have shown that attentional facilitation lingers at the retinotopic coordinates of a previously attended position after an eye movement. These results are intriguing, because the retinotopic location becomes behaviorally irrelevant once the eyes have moved. Critically, in these studies participants were asked to maintain attention on a blank location of the screen. In the present study, we examined whether the continuing presence of a visual object at the cued location could affect the allocation of attention across eye movements. We used a trans-saccadic cueing paradigm in which the relevant positions could be defined or not by visual objects (simple square outlines). We find an attentional benefit at the spatiotopic location of the cue only when the object (the placeholder) has been continuously present at that location. We conclude that the presence of an object at the attended location is a critical factor for the maintenance of spatial constancy of attention across eye movements, a finding that helps to reconcile previous conflicting results

Re-sampling strategy to improve the estimation of number of null hypotheses in FDR control under strong correlation structures

<p>Abstract</p> <p>Background</p> <p>When conducting multiple hypothesis tests, it is important to control the number of false positives, or the False Discovery Rate (FDR). However, there is a tradeoff between controlling FDR and maximizing power. Several methods have been proposed, such as the q-value method, to estimate the proportion of true null hypothesis among the tested hypotheses, and use this estimation in the control of FDR. These methods usually depend on the assumption that the test statistics are independent (or only weakly correlated). However, many types of data, for example microarray data, often contain large scale correlation structures. Our objective was to develop methods to control the FDR while maintaining a greater level of power in highly correlated datasets by improving the estimation of the proportion of null hypotheses.</p> <p>Results</p> <p>We showed that when strong correlation exists among the data, which is common in microarray datasets, the estimation of the proportion of null hypotheses could be highly variable resulting in a high level of variation in the FDR. Therefore, we developed a re-sampling strategy to reduce the variation by breaking the correlations between gene expression values, then using a conservative strategy of selecting the upper quartile of the re-sampling estimations to obtain a strong control of FDR.</p> <p>Conclusion</p> <p>With simulation studies and perturbations on actual microarray datasets, our method, compared to competing methods such as q-value, generated slightly biased estimates on the proportion of null hypotheses but with lower mean square errors. When selecting genes with controlling the same FDR level, our methods have on average a significantly lower false discovery rate in exchange for a minor reduction in the power.</p

Effects of dependence in high-dimensional multiple testing problems

<p>Abstract</p> <p>Background</p> <p>We consider effects of dependence among variables of high-dimensional data in multiple hypothesis testing problems, in particular the False Discovery Rate (FDR) control procedures. Recent simulation studies consider only simple correlation structures among variables, which is hardly inspired by real data features. Our aim is to systematically study effects of several network features like sparsity and correlation strength by imposing dependence structures among variables using random correlation matrices.</p> <p>Results</p> <p>We study the robustness against dependence of several FDR procedures that are popular in microarray studies, such as Benjamin-Hochberg FDR, Storey's q-value, SAM and resampling based FDR procedures. False Non-discovery Rates and estimates of the number of null hypotheses are computed from those methods and compared. Our simulation study shows that methods such as SAM and the q-value do not adequately control the FDR to the level claimed under dependence conditions. On the other hand, the adaptive Benjamini-Hochberg procedure seems to be most robust while remaining conservative. Finally, the estimates of the number of true null hypotheses under various dependence conditions are variable.</p> <p>Conclusion</p> <p>We discuss a new method for efficient guided simulation of dependent data, which satisfy imposed network constraints as conditional independence structures. Our simulation set-up allows for a structural study of the effect of dependencies on multiple testing criterions and is useful for testing a potentially new method on <it>π</it>₀or FDR estimation in a dependency context.</p

Scans for signatures of selection in Russian cattle breed genomes reveal new candidate genes for environmental adaptation and acclimation

Domestication and selective breeding has resulted in over 1000 extant cattle breeds. Many of these breeds do not excel in important traits but are adapted to local environments. These adaptations are a valuable source of genetic material for efforts to improve commercial breeds. As a step toward this goal we identified candidate regions to be under selection in genomes of nine Russian native cattle breeds adapted to survive in harsh climates. After comparing our data to other breeds of European and Asian origins we found known and novel candidate genes that could potentially be related to domestication, economically important traits and environmental adaptations in cattle. The Russian cattle breed genomes contained regions under putative selection with genes that may be related to adaptations to harsh environments (e.g., AQP5, RAD50, and RETREG1). We found genomic signatures of selective sweeps near key genes related to economically important traits, such as the milk production (e.g., DGAT1, ABCG2), growth (e.g., XKR4), and reproduction (e.g., CSF2). Our data point to candidate genes which should be included in future studies attempting to identify genes to improve the extant breeds and facilitate generation of commercial breeds that fit better into the environments of Russia and other countries with similar climates

The Mitochondrial Genome of the Legume Vigna radiata and the Analysis of Recombination across Short Mitochondrial Repeats

The mitochondrial genomes of seed plants are exceptionally fluid in size, structure, and sequence content, with the accumulation and activity of repetitive sequences underlying much of this variation. We report the first fully sequenced mitochondrial genome of a legume, Vigna radiata (mung bean), and show that despite its unexceptional size (401,262 nt), the genome is unusually depauperate in repetitive DNA and "promiscuous" sequences from the chloroplast and nuclear genomes. Although Vigna lacks the large, recombinationally active repeats typical of most other seed plants, a PCR survey of its modest repertoire of short (38–297 nt) repeats nevertheless revealed evidence for recombination across all of them. A set of novel control assays showed, however, that these results could instead reflect, in part or entirely, artifacts of PCR-mediated recombination. Consequently, we recommend that other methods, especially high-depth genome sequencing, be used instead of PCR to infer patterns of plant mitochondrial recombination. The average-sized but repeat- and feature-poor mitochondrial genome of Vigna makes it ever more difficult to generalize about the factors shaping the size and sequence content of plant mitochondrial genomes

Thyrotropin-releasing hormone (TRH) promotes wound re-epithelialisation in frog and human skin

There remains a critical need for new therapeutics that promote wound healing in patients suffering from chronic skin wounds. This is, in part, due to a shortage of simple, physiologically and clinically relevant test systems for investigating candidate agents. The skin of amphibians possesses a remarkable regenerative capacity, which remains insufficiently explored for clinical purposes. Combining comparative biology with a translational medicine approach, we report the development and application of a simple ex vivo frog (Xenopus tropicalis) skin organ culture system that permits exploration of the effects of amphibian skin-derived agents on re-epithelialisation in both frog and human skin. Using this amphibian model, we identify thyrotropin-releasing hormone (TRH) as a novel stimulant of epidermal regeneration. Moving to a complementary human ex vivo wounded skin assay, we demonstrate that the effects of TRH are conserved across the amphibian-mammalian divide: TRH stimulates wound closure and formation of neo-epidermis in organ-cultured human skin, accompanied by increased keratinocyte proliferation and wound healing-associated differentiation (cytokeratin 6 expression). Thus, TRH represents a novel, clinically relevant neuroendocrine wound repair promoter that deserves further exploration. These complementary frog and human skin ex vivo assays encourage a comparative biology approach in future wound healing research so as to facilitate the rapid identification and preclinical testing of novel, evolutionarily conserved, and clinically relevant wound healing promoters

Arabidopsis MKS1 Is Involved in Basal Immunity and Requires an Intact N-terminal Domain for Proper Function

Innate immune signaling pathways in animals and plants are regulated by mitogen-activated protein kinase (MAPK) cascades. MAP kinase 4 (MPK4) functions downstream of innate immune receptors via a nuclear substrate MKS1 to regulate the activity of the WRKY33 transcription factor, which in turn controls the production of anti-microbial phytoalexins.We investigate the role of MKS1 in basal resistance and the importance of its N- and C-terminal domains for MKS1 function. We used the information that mks1 loss-of-function partially suppresses the mpk4 loss-of-function phenotype, and that transgenic expression of functional MKS1 in mpk4/mks1 double mutants reverted the mpk4 dwarf phenotype. Transformation of mks1/mpk4 with mutant versions of MKS1 constructs showed that a single amino acid substitution in a putative MAP kinase docking domain, MKS1-L32A, or a truncated MKS1 version unable to interact with WRKY33, were deficient in reverting the double mutant to the mpk4 phenotype. These results demonstrate functional requirement in MKS1 for the interaction with MPK4 and WRKY33. In addition, nuclear localization of MKS1 was shown to depend on an intact N-terminal domain. Furthermore, loss-of-function mks1 mutants exhibited increased susceptibility to strains of Pseudomonas syringae and Hyaloperonospora arabidopsidis, indicating that MKS1 plays a role in basal defense responses.Taken together, our results indicate that MKS1 function and subcellular location requires an intact N-terminus important for both MPK4 and WRKY33 interactions